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An alternative parametrization of particle shower longitudinal profile is presented. The accuracy 
of obtained shower profile description is about 2 — 3% for the 0 — 1500 g/cm 2 atmosphere slant 
depths and primary H, He,.. . Fe nuclei in 1 PeV-10 EeV energy range. It is shown that the shape of 
shower profile depends only on the nucleon energy, whereas the maximum shower size also depends 
on the energy of parental nucleus. Results are based on the CORSIKA simulated shower profiles and 
are presented in comparison with Gaisser-Hillas parametrization. 

PACS numbers: 96.50.S-, 96.50.sd, 96.50.sb 


I. INTRODUCTION 

The parametrization of longitudinal profiles for par¬ 
ticle showers produced by the primary nuclei in atmo¬ 
sphere, is an essential tool for the primary nuclei iden¬ 
tification and the evaluation of primary energy. Experi¬ 
ments that sample the shower longitudinal development 
using Cherenkov light images [l] or air fluorescence Hi 
from different traversed atmosphere depths, extract the 
position of shower maximum, which is sensitive to the 
incident primary nucleus. The integral of shower profile 
strongly correlates with primary energy Q. 

The shower longitudinal profile is a dependence of the 
shower particle number (N) on a given traversed atmo¬ 
spheric depth, T. The parametrization of shower pro¬ 
file commonly used in cosmic-ray experiments is Gaisser- 
Hillas formula [fij: 

N(X) = N max ( -2—) exp (. X max - X) , (1) 

where X = (T - X 0 )/ A and X max = ( T max - X 0 )/A. 

The maximum number of shower particles N max at the 
traversed atmosphere depth T max along with Ao and A 
in expression (1) are free parameters that depend on the 
primary nucleus and energy. 

The standard primary nuclei composition consists of 
the first 28 nuclei of the periodic table with mass (nu¬ 
cleon) numbers A = 1,... 56 usually divided into four-six 
groups (species) H, He, CNO-like, Si-like, Fe-like. The 
large number of nuclei species (more than four) increases 
the uncertainties of the inverse problem (E and A recon¬ 
struction) falsely improving the agreement of experiment 
with the theory [6j. 

The primary energy region responsible for particle 
shower detection at the observation level begins at about 
E > 1 PeV and ends at GZK cutoff energies Q- 

The efficiency of four-parametric parametrization (1) 
is in its applicability to a wide range of energies and pri¬ 
mary nuclei. However, the observed correlations between 
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parameters result in a loss of the physical meaning of Xq 
and A [?] and reduce the range of effective atmosphere 
depths for Eq. (1). 


II. PARAMETRIZATION 

Here, an alternative parametrization N(T , E, e) for 
particle shower longitudinal profile is proposed using 
three non-correlating parameters that depend on the pri¬ 
mary particle energy and nucleon energy, e: 

N{x) = A max exp(—-1 ) > ( 2 ) 

where 

8(x) = a — /3(tanha;) 3 (3) 

is the profile shape function of variable 

T 


The shower maximum position, T max (e), and shape 
function, 8(x, e), turned out to be dependent on the pri¬ 
mary particle energy per nucleon, 

£ = \ , (PeV/n). 

The maximum number of shower particles, N max (E, e) 
is factored into the primary energy and a function of 
nucleon energy only. The corresponding approximations 
for the parameters of shower longitudinal profile (2,3) 
are: 

a = 0.707 + 0.209e'°- 084 , /3 = yV2.59 , (4a) 

Tmax = 433.5 + 38.9(ln(eAi? e ))°' 857 , (g/cm 2 ), (4b) 

N m ax = 0.653(£7/lGeV)(l - e- 2 ' 5£ ° 12 ) , (4c) 

where Ap e = 56 and e is in the units of PeV/n. The 
goodness-of-fit tests for (4a-c) were y 2 < 1 at negligible 
correlations between the a, T ma x and N max parameters. 
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FIG. 1: Average longitudinal shower profiles for 13 traversed 
atmosphere depths produced by H, He, O and Fe primary 
nuclei with 6 energies (from 1 PeV to 10 EeV). The symbols 
are CORSIKA shower simulated data (training sample). The 
dashed lines are results from the 4-parametric approximation 
(1). The solid lines are parametrizations (2-4) computed for 
corresponding primary nuclei and energies. 


III. SHOWER PROFILES 

The values of free parameters in expressions (1) and 
(2) were obtained from simulated shower profiles (train¬ 
ing sample) using CORSIKA 0 (SIBYLL 1 ) code for 
four primary nuclei A = 1,4,16,56 at six energies 
E = 1,10,100, 500, 2500,10 4 PeV. Shower profiles were 
studied for 10 atmosphere depths T = 100, 200,... 1000 
g/cm 2 at two zenith angles, cosd = 0.7 and 1. The 
shower particle energy threshold was E e > 1 MeV. Sim¬ 
ulation statistics were provided for less than 2-3% statis¬ 
tical errors in the whole measurement range. 

The averaged shower profiles were approximated by 
expressions (1) and (2) using 13 reference depths. The 
results are presented in Fig. [1] It is seen that the 
parametrization (1) (dashed lines) underestimates the 
shower sizes at large atmosphere depths. 

The parametrization errors of expressions (1) and (2) 
and corresponding \ 2 d o f are presented in Fig. [2] for 
different primary energies and nuclei. The upper and 
middle panels show the errors of the four-parametric ap¬ 
proximations of CORSIKA simulated shower profiles using 
the N max , T max , Xo and A parameters of expression (1) 
and the N max ,T max , a and /3 of expressions (2,3). The 
lower panel of Fig. [2] shows the errors of shower profiles 
N(T,E,A) from expressions (2-4). 

The normalized simulated (symbols) and parametrized 
(lines) shower profiles are presented in Fig. [3] It is 
seen that the parametrization (2) effectively describe 
the shower profiles in the regions of both the maximum 
(x ~ 1, inset figure, solid line) and asymptotic depths 


FIG. 2: Parameterization errors ( Nf it — N cors ika) / AN corsi ka 
for the 4-parametric expressions (1) and (2,3) are shown in the 
upper and middle panels respectively. The lower panel cor¬ 
responds to the errors of the shower profiles N(T, A, E) from 
expressions (2-4) for different primary energies and nuclei. 



FIG. 3: Normalized particle shower profiles. The symbols are 
CORSIKA simulated data (training sample). The lines are 
the results of parametrization (2). Inset figure is a zoom of 
selected rectangular region for 0.8 < x < 1.2. The dashed 
lines correspond to the parametrization (1). 


{x ~ 3). Equation (1) is systematically biased about 
—2% (inset figure, dashed line) at x ~ 1. 


IV. PARAMETERS 

The study of T max (e) and N max (e) dependence on nu¬ 
cleon energy (e) are presented in the upper and lower 
panel of Fig. [I] respectively. The approximations of 
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FIG. 4: Parameters Tmax (upper panel) and N max /E in units 
of particle/GeV (lower panel) derived from expressions (1) 
and (2) for the different boundaries of traversed depths. The 
lines correspond to the expressions (4b) and (4c) for T max and 
Nmax respectively. 

shower profiles using parametrizations (1) and (2) were 
trailed for different lower ( Ti ow ) and upper (T up ) limits 
of traversed atmosphere depth. 

The estimated values for T max (Fig. [4] upper panel) 
were unbiased for all trails. The line in Fig. [4] corre¬ 
sponds to the expression (4b). The asterisk and cross 
symbols in Fig. [5] are correspondingly renormalized COR- 
SIKA simulated data from [lOJ. 

Estimations of N rnax (Fig. [2 lower panel) using ex¬ 
pression (1) for approximations of shower profile turned 
out to be dependent on boundary conditions for atmo¬ 
sphere depth (hollow and bold star symbols), whereas 
the expression (2) remained practically unbiased (hollow 
and bold circle symbols) for different boundaries. 

The shower profile shape functions <5(e|T) and <5(cc|er) 
are presented in Fig. [5j where the symbols (left panel) 
are the data extracted from CORSIKA simulated training 
sample. The solid lines in both panels correspond to the 
expressions (3,4a). The dashed lines in the right panel of 
Fig. [5] are the 0.5% accuracy logarithmic simplifications 
of shape function (3), 

a — b\ ni, if 0.07 < x ^ 1; 

a — b\nx/{\ + 1.59 In a;), if x ^ 1, 

(5) 

where 

a = 0.215+ 0.145e~ 0 084 , 

' /M 

b = 0.086 + O.Olle -0 084 , 

at x 2 /450 = 0.7. The approximation (5) provides the 



e=E/A, (PeV/n) Atm. depth, x=T/T mox 

FIG. 5: Shower profile shape function (expression (3)) de¬ 
pending on the nucleon energy (left panel) and normalized 
atmosphere depth (right panel). The solid lines correspond 
to the expressions (3,4a). The dashed lines in the right panel 
are the logarithmic simplifications of the shape function ac¬ 
cording to (5,6). 

analytic solution of the inverse function Nq (x) for x < 1 
and N^ 1 (x) for x > 1. 


V. VERIFICATION 

The verification of the universality of approximations 
(1) and (2) was performed by extrapolating the shower 
profiles from 100 — 1428 g/cm 2 interval to the T = 10 
g/cm 2 observation level, corresponding to the earlier 
stage of shower development. The results are presented 
in Fig. [G] The symbols at T = 10 g/cm 2 in Fig. [6] are 
the corresponding data from CORSIKA simulated control 
sample, whereas the symbols at T = 100 g/cm 2 are the 
representatives of the training sample (Section III). 

It is seen that parametrization (1) being trained in the 
100 — 1428 g/cm 2 depth interval can not be extrapolated 
to the region less than about 50 g/cm 2 (lines, left panel), 
whereas the parametrization (2) works correctly up to 
the beginning of the atmosphere (lines, right panel). 

The verification of the shower profiles (2-4) by the con¬ 
trol samples of different nuclei and energies are shown in 
Fig. [3 The shower profile for primary Fe nucleus with 
energy E = 500 PeV and corresponding e ~ 8.93 PeV/n 
from training sample (Section III) are compared to the 
control sample of shower profiles produced by primary H, 
He, C, O and Si nuclei with the same energy per nucleon 
(symbols). The lines in Fig. [3 are the corresponding con¬ 
gruent predictions from the parametrizations (2-4). 

The results in Fig. [3 confirm the e-dependence of the 
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with the accuracies of about 2 — 3% in the whole mea¬ 
surement range. 


VI. INTEGRAL 

The right hand side of parametrization (2) at the cor¬ 
responding normalization can be considered as a prob¬ 
ability density function and be used for primary energy 
evaluation [2|, 0] - Unfortunately this function was missed 
by mathematicians and by using numerical technique the 
required normalization 

f(x, e)da; ~ 1 ± 10 -4 (7) 

was provided for probability density function 




FIG. 6: Extrapolations of the parametrization (1) (left panel) 
and (2) (right panel) to the earliest stage of shower develop¬ 
ment for different primary nuclei and energies (lines). The 
symbols are the CORSIKA simulated data. 
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with additional parameter 

do = 0.226+ 0.148e-° 092 


The goodness-of-fit test for do(e) was % 2 = 0.01 in the 
1CT 2 < e < 10 4 (PeV/nucleon) interval and the upper 
limit of integral (7), x max = 3. 

It is interesting to note the relation between parame¬ 
ters do and shape function d(x) from expression (3): 


d 0 (e) ~ 


d(a:)dx ± 1% . 


(9) 


The statistical parameters, the average ( x) and standard 
deviation (a x ) of distribution (8), are well approximated 
(0.1% errors) by the following expressions that depend 
on the nucleon energy: 


x = 1.036+ 0.094£-° 12 


at x 2 = 0.1, and 

a x = 0.226+ 0.176£-°° 92 


at x 2 = 1-1- 


FIG. 7: Control samples of shower profiles (symbols) pro¬ 
duced by the different primary nuclei with the same nucleon 
energy e ~ 8.93 PeV/n. The lines are the predictions from 
(2-4). 


shower longitudinal profile shape (expressions 4a,b). The 
shower profile amplitude (N max ) also depends linearly on 
the primary energy, E (expression 4c). 

The good agreement in Fig. [7] between predictions 
(lines) and simulated data indicates the correctness of 
expressions (2-4) for shower profile description at least 


VII. FLUCTUATIONS 

The main source of shower profile fluctuations is the 
depth of the first interaction of primary particles in the 
atmosphere The exponentially distributed uncer¬ 

tainty of the first interaction point results in the corre¬ 
sponding fluctuations of the shower profile (2) depend¬ 
ing on the rate of change (dlV/da:) of the profile with 
respect to the depth, x. Thus, the fluctuations should 
be maximal at the beginning of shower development 
(x ~ 0, Fig. [3]), and minimal in the region of shower 
maximum, x = 1. The dependence of the interaction 
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FIG. 8: Normalized standard deviations ( <jn/N ) of shower 
particles for different primary nuclei and primary energies 
(symbols). The lines are represent the parametrization (10) 
for energies of 1 PeV (dotted lines) and 10 EeV (solid lines). 
The dashed lines describe the fluctuations for intermediate 
energies. The inset panel zooms in on the region of minimal 
fluctuations at x ~ 1. 

length, A (A,E), on the primary particle also results in 
the mass (A) and energy ( E ) dependences of shower pro¬ 
file fluctuations. 

The statistical measure of fluctuations is the standard 
deviation of shower particle number, cnv- The corre¬ 
sponding values of <jn{x, A, E)/N(x) obtained from the 
shower simulated dataset (Section III) are presented in 
Fig. [5] (symbols). The inset panel shows the region of 
minimal fluctuations in detail. The lines in Fig. [8] corre¬ 
spond to the parametrizations 

&n_ ~ j a-! - a 2 lnx, if x ^ 1; , , 

N |ai + a 3 (ln' ? x)/x, if x ^ 1, 

where 

ai = 0.165A _0 ' 32 .E _013 , 
a 2 = 0.68A-°- 185 Fr o o ° 9 , 
a 3 = 3 . 77 ^- 0 - 386 ^- 0.035 _ 

r] = 2.67A-°- 080 E-° 027 

at x 2 /470 ~ 1.7. 

VIII. SUMMARY 

The standard inverse problem of cosmic ray physics in 
PeV-EeV energy region is the identification of a primary 


nucleus (or elemental composition) and the estimation of 
its energy (or energy spectrum) by the detected shower 
response at the observation level in the frames of a given 
interaction model. The efficiencies of primary particle 
and primary energy estimators depend on both the accu¬ 
racy (Section III) and universality (Section V) of shower 
longitudinal profile description. 

Historically, the conventional shower longitudinal pro¬ 
files were proposed in 1960 (Greisen function) , 1977 
(Gaisser-Hillas function) 15| and 2001 (Gaussian-In-Age 
approach) [Ti|. The efficiencies and accuracies of listed 
profile parametrizations are compared in Refs. il. IT tL [l5l] 
in detail. 

The last Gaussian-In-Age approach Q reduced the 
number of parameters to 3, and decreased the inter¬ 
correlations between parameters of profile function 
in return for the narrow range of applicability in the 
vicinity of shower maximum: 0.75 < s < 1.25 [14j . where 
s = 3/(1 + 2/x) is the shower age parameters. 

The alternative shower longitudinal profile descrip¬ 
tion (expressions 2-4), as opposed to the parametriza¬ 
tions s i, m , represents the first complete formula 
for shower profile, N(T, A, E), depending on the atmo¬ 
sphere depth ( T ), primary nucleus ( A ) and primary en¬ 
ergy E. Expressions (2-4) provide the accuracies of about 
2 — 3% for the region of 0 < T < 1450 g/cm 2 , A < 56, 
1 PeV < E < 10 EeV. The results are obtained in the 
frames of SIBYLL 0] interaction model (Section III). 

The position of the shower maximum, T max (£) from 
expression (4b) and profile shape function, S(x,£ ) from 
expression (4a), depend only on the primary nucleon en¬ 
ergy £ = E/A , which is in agreement with the prediction 
of superposition model 0. 

The amplitude of the profile N max (E,£) from expres¬ 
sion (4c) depends on both the primary energy (E) and 
nucleon energy (e). 

The intercorrelations between the N max (E 1 £), 
Tmax(s) and <5( x,e) shower profile parameters are 
negligible. 

The profile shape function, 5(x, e), from (3) has the 
simple logarithmic representation (5) which provides an 
analytic solution for the corresponding inverse profile 
function, which can be used in the Constant-Intensity- 
Cut method Q. 

The fluctuations of particle shower longitudinal profile, 
<t/v/A 7, from parametrization (10) depend on the energy 
(E) and mass number ( A ) of the primary nuclei. 
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